Equation of state of a polarized Fermi gas in the Bose-Einstein Condensate limit 



F. Alzetto and X. Leyronas 
Laboratoire de Physique Statistique, Ecole Normale Superieure, UPMC Paris 06, 
Universite Paris Diderot, CNRS, 24 rue Lhomond, 75005 Paris, France. 

We present a theoretical study of the BEC-BCS crossover in the Bose-Einstein-Condensate regime 
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as the one-particle gap and the energy of an "impurity" immersed in a Fermi sea. For the chemical 
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I. INTRODUCTION 

The problem of a Bose-Fermi mixture in the Bose- 
Einstein-Condensate (BEC) limit of the BEC-BCS 
(Bardeen-Cooper-Schrieffer) crossover has been studied 
experimentally recently pQ. Previous theoretical works 
studied the zero temperature phase diagram using a 
mean field approach [2], or Monte Carlo simulations [3]. 
The BEC limit has been investigated theoretically in 
Refs.[4],[5] or using a functional integral formalism [6]. 
In Ref. [|] , the authors used a Bogoliubov-de Gennes ap- 
proach in order to derive in an elegant way a Gross- 
Pitaevskii equation for the condensate wave function. 
They found a dimer-dimer scattering length and an atom- 
dimer scattering length equal to 2 times the atom-atom 
scattering length a and 8/3 a respectively. Then they 
gave some reasonable diagrammatic arguments in order 
to recover the correct scattering lengths. As far as the 
equation of state of the system is concerned, these calcu- 
lations stopped at the mean-field level (where the chem- 
ical potentials depend linearly on the densities). On the 
other hand, researchers had studied previously [7], [8] a 
mixture of point-like bosons and fermions at zero tem- 
perature in the dilute limit. There the equation of state 
was determined including a lowest order mean-field cor- 
rection. This is the main result of the present work to re- 
cover this result for the problem of a mixture of fermions 
and composite bosons relevant to experiment of Ref. pQ . 

In this work we determine the equation of state at zero 
temperature of a mixture of fermions of two species. We 
consider a dilute situation where the densities of both 
species n\ and ni (hereafter denoted density) are small 
but are of the same order. We thus have density x a 3 <C 
1, and we expand the chemical potentials of fermions up 
to order (density) 4 / 3 . 

For a reader mainly interested in the results, and not 
so much in technical details, the main results of our work 
(the equations of state) are shown in Eqs 51|54 and 55 



In Sec|IIJ we recall the general formalism of fermionic 
superfluids (anomalous Green's functions, anomalous 
self-energy, and the so-called generalized self-energies). 



In Sec |III[ we describe the general framework of the the- 
ory, namely how we expand the number and gap equa- 
tions in powers of the atomic density. In Sec |IV[ the low- 
est order density expansion is deduced, while the deriva- 
tion of the next order correction is shown in SecJV] The 
one-particle gap and the problem of a single atom im- 
mersed in the Fermi Sea of atoms of the other species 
are determined in Sec |VI[ In Sec |VII[ we conclude and 
comment on the relevance of our calculation for recent 
experiments of Ref. [T] . 



II. GENERAL FORMALISM 

In a fermionic superfluid, the usual field-theoretical 
methods must be modified in order to take into account 
the off-diagonal order [S]. In practice, one has to consider 
"anomalous" Green's functions F(p) and F + (p) (which 
are physically the propagators of condensed pairs). The 
Dyson equations for the fermionic superfluid are 

Gi( P ) = G?(p) + GSOOExOOGiOO + G° 1 ( P )A(p)F + (p) 

(1) 

for the propagator of particles of type 1 (and similarly 
for particles of type 2) and 

F+(p) = Gl(-p)^(-p)F+{p)-Gl{-p)A*{p)G 1 {p) 

(2) 

for the anomalous Green's function F + (p). These two 
equations have a simple diagrammatical meaning (see 
Figjl]): the self-energy £i(p) is the sum of all diagrams 
with one incoming (p, 1) line (a G\(p) Green's function) 
and one outgoing (p, 1) line, which can not be separated 
into two disconnected parts by cutting a (p, 1) line or a 
(— p, 2) line (a G^i—p) Green's function). The anoma- 
lous self-energy A(p) is the sum of all diagrams with one 
incoming (p, 1) line and one incoming (— p, 2) line, which 
can not be separated into two disconnected parts by cut- 
ting a (p, 1) line or a (— p, 2) line. The anomalous Green's 
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function F + (p) starts diagrammatically with one outgo- 
ing (— p, 2) line and finishes with one outgoing (p, 1) line. 
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FIG. 1: The Dyson equations Eq[TJ and [2] in diagrammatic 
language. The numbers 1 and 2 specify the atomic species. 



In Ref.[TT], the so-called generalized self-energy Si(p) 
was introduced. It is defined by the equation 



Gi(p) = 



(G2(p))-i-Si(p) 



(3) 



which can be rewritten G\(p) = G\(p) + 

G!(p)^i(p)G?(p) + G^SiiP^MSMG^p) + •••• 
Diagrammatically, this means that Si(p) is the sum 
of all diagrams with one incoming (p, 1) line and one 
outgoing (p, 1) line which cannot be disconnected into 
two parts by cutting a (p, 1) line. It is natural to use 
the generalized self energy, since it is the only thing one 
needs in order to calculate the one particle Green's func- 
tion Gi(p) (whereas one needs both and A(p)) and 
therefore the density of fermions of species 1. Moreover, 
as it was shown in Ref. , the generalized self-energy is 
directly related to the dimer-atom T-matrix (hereafter 
denoted T 3 , as we will see in Section IV), whereas £i(f>) 
is not. This makes the use of the generalized self-energy 
convenient to get more compact expressions. 



III. NUMBER AND GAP EQUATIONS 

The system we study is an ensemble of two distinct 
"species" of fermions. There are n\ (respectively nz) 
fermions per unit volume of type 1 and mass mi (resp. 
type 2 and mass m-z). These two species can be either two 
internal states (for instance two hyperfine states of 6 Li) 
and in this case m\ = 7712 or two different sort of fermions 
(a 6 Li- 40 K mixture for instance). Due to the extremely 
low temperature and Pauli exclusion principle there is no 
intra-species interaction. The short range inter-species 
interaction is therefore described by the scattering length 
a (for a review, see Ref. [TO]). In the limit we consider, 
we have a > and there is a two-body bound state (a 
molecule or dimer) of energy E\, = —l/(2m r a 2 ) where 
m r = m\mij(rn\ + 777,2) is the particle 1-particle 2 re- 
duced mass and we will take in this work K = 1. As 
pointed out by Pieri and Strinati @], the limit of dilute 
condensed molecules with extra fermions of species I cor- 
responds to a positive chemical potential fii and a nega- 
tive and large chemical potential (of order Ej,) fj, 2 for the 
fermions 2. Therefore this is the limit we will consider 



too and we take 



and 



(H > 
M2 < 

Hi + (12 < 



(4) 
(5) 

(6) 



The method is similar to the one used in Ref. [TTJ and [T2] , 
and is common to most of the studies of the BCS-BEC 
crossover |I0j : we calculate the particles number ni and 
ni as a function of the chemical potentials /ii and \xi 
and a superfluid order parameter A; we write down also 
a generalized "gap equation" for A as a function of /ii 
and [12 ■ If we eliminate A from the equations, we get the 
equations of state /zi (711,712) and ^(ni, n?*)- The results 
for these equations of state are shown in Eqs |51|54| and 
|55[ which are the main results of the present work. 



A. Number equations 

In order to get the density of particle 1, we have to 
sum on k the propagator Gi(k) 



ni 



(7) 



where we use the notation J^k = — / rjL) 3 ^ S n ( m an 
the paper, we use the T ^ formalism and then take 
the limit T ->• 0) . We have G\(k) = -, x - — 

with k = {k, iio n }, fii > and uj n — (2n + 1)ttT (with 
n integer). We have of course a similar result for the 
Green's function of atoms 2, with a chemical potential 
/12 < 0. Following (TTJ, the full Green's function Gi(k) is 
related to the generalized self-energy Si(k) via a Dyson 
equation Gi(fc) = G?(fc) + G° 1 (k)Si(k)G 1 (k). Si(k) does 
not contain diagrams with a G\(k) propagator, but it 
does contains diagrams with G^{— k) propagators. If wc 
iterate the Dyson equation, we get the approximation 



Gi(k) ^G°i(k) + (G°i(k)) 2 Si(k) 



(8) 



As we explain below, this approximation is good enough 
at the order we consider. It must be emphasized that the 
generalized self-energy contains anomalous contributions, 
as we will see. Using the approximation EqjHJ we get 
from the first term the contribution to the density (after 
taking the T — > limit) 



1 6tt 2 

where the wave-vector fco is defined via 
fc = ^/2mini 



If we define Srii 



ii — rii, we get 



8m « -£(G?(fc)) a £i(fc) 



(9) 
(10) 

(11) 
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We can proceed in the same way for the density of par- 
ticles 2. In this case we have n\ = 0, since \ii < 0. We 
therefore have 



»2 



-£(G§(£0) 2 S 2 (fc) 



(12) 



where S2(k) is the generalized self-energy of the particles 
2. 



for TZe[iuj u ) < (in many instances in this work, we find 
convenient to use the somewhat loose notation g[iuj v ) 
for the analytical continuation of a function g for iu> v 
spanning the whole complex plane - not only the discrete 
Matsubara frequencies values). The final gap equation 
reads 



m, 
~2n 



t fa' 1 - y^mM) = £>' (fc)/A) - i?(0)(18) 



B. Gap equation 

We follow Ref.[TT] in order to get our gap equation. 
Looking at the diagrammatic expansion of the anomalous 
Green's function F[k), we can in general write : 

F(k) = G?(fc)G°(-fc) ^^(k-kOFfo)^ 



+F'[k) 



(13) 



The first term in the r.h.s. of Eqfl3| contains all the dia- 
grams where the two incoming (As, 1) and (—k, 2) fermions 
lines first interact (interaction V(q}). Therefore F'(k) 
is the sum of all the diagrams where the two incoming 
fermions (As, 1) and (—As, 2) do not first interact. The term 
between parentheses in Eq{l3] is the usual BCS gap [TT] 
Mk) 

£r(k) = J2 V ^-^i)F(ki) (14) 

fci 

We can eliminate the interaction V in favor of the 2-body 
T-matrix [llj at zero momentum T 2 (0). 



[^(O)]- 1 = $>'(As)/A 



(15) 



where A is the zero wave vector value of Si (k) . As usual, 
Ti is obtained by summing the ladder series and one finds 
easily 

T 2 [K)~ 1 = g^-UiK) 



2tt 



— 2m r (£f2„ + 2/i — 



K 2 

2 mjj ' 



-R[K) 



(16) 



where K = {K,iu! u }, uj v — 2ttvT [v integer), mu = 
mi + to 2 is the mass of a dimer and 2/x = jjL\ + M2 is 
the total chemical potential. R{K) is due to the positive 
chemical potential fii : 



R(K) = 



d 3 k 



1 



k<k Q (2tt) 3 f iQu + 2^ 



_k? (K-k) 2 

2 mi 2 m,2 



(17) 



This is due physically to the presence of the Fermi Sea of 
the fermions in excess. We note that T2(K), is analytical 



with R(0) = - f, ,fX t u3 1 ^ 

V I Jk<k„ (2tt) 3 fj£-+2\n\\ 

We define the bosonic chemical potential \xb through 
the equation 



Hi + H2 = 2/x 



\E b \+HB 



(19) 



and we can expand the l.h.s. of Eq |18| in power of hb- In 
this way we get 



Mb 



C£F>(k)/A)-R(0) 



O(hI) 



(20) 



where we have defined the residue of T 2 at the molecular 
pole 



2tt 



(21) 



In the following, we will expand F'(k), Si(k) and 52 (As) 
in power of A and Aso. We will see this corresponds to an 
expansion in power of the density of the condensed dimer 
and of the fermions in excess respectively. 



IV. DENSITY EXPANSION : LOWEST ORDER 

A. Number equations 

We start with the calculation of Si (As). Diagrams with- 
out anomalous self-energy are diagrams in the normal 
phase. In this case, one has to close a fermionic loop 
for the species 2. Going back to time representation, 
one sees that such contributions are zero, since the bare 
propagators of particles 2 are retarded, as is explained in 
details in [TT]. Physically, this comes from the fact that 
at T = 0, there are no "free" [i.e. not bound to a con- 
densed molecule) fermion of species 2 which can interact 
with a fermion of species 1. The diagrams of lowest order 
in A (that is of lowest order in y/n^) are therefore those 
containing one A and one A*. They are given in terms 
of the 3-particle T-matrix Tg 1 describing the scattering of 
a dimer with fermion of species 1 (see Figj2ja)) 



S[ a \k) = \A\ 2 Ti(k,k;k) 



(22) 



This diagram has a very simple physical meaning. 
An incoming fermion k, 1 scatters with a condensed 
dimer (factor A* for a dimer "out of the conden- 
sate), the scattering of the dimer with the fermion 
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is described by the dimer-atom 1 T-matrix Tg 1 . The 
right of the graph represents the outgoing fermion fc, 1 
and the dimer which "goes back" to the condensate 
(factor A). The contribution to the density of par- 



ticles 1 is 8n\ 



(a) _ 



(G°(fc)) 2 |A| 2 T 3 HM;fc). 



The frequency summation is written as usual as 
an integral on a contour encircling the imaginary 
axis counterclockwise : TJ2 iu JG^(k)) 2 T^(k, k; k) = 

~i^f( z )(^ T i(^ k ^) with/(z) = l/( e ^ + l)is 
the Fermi Dirac distribution and £i ( 2 = k 2 /(2mi i 2) — /ii.2 
and k — {k, z}. We use next the analytical properties of 
T3 as a function of z. Following an argument given in 
[IT] , we can split T 3 in a first part l/(z + £ 2 ) which has 
a pole at —£2 < plus a function which is analytical for 
IZe(z) < 0. When we deform the contour in order to en- 
close the singularities, we see that the only singularities 
in IZe(z) < are — £ 2 and £1 if £1 < 0. Due to the ex- 
ponential decay of the Fermi-Dirac distribution f(z) for 
IZe(z) < when the temperature goes to 0, the contour 
encircling the domain IZe(z) > vanishes. We therefore 
find at zero temperature 



Tj2G°i(kfTi(k,k;k) 



(C1 + 6) 2 



dz 



+/(6)^T 3 1 (fc,fc;fc) (23) 



where the 4- wave vector k = {k, £1} is "on the shell". In 
order to compute the density the next step is to integrate 
the three terms on the r.h.s. of Eq(23]on the wave vector 
k. We use the zero temperature limits ^jf 1 ^ = — <K£i), 
= O(fco ~ |k|) is the step function) and 



/(— £2) = 1- The first term is simply f-M^ 



+2|H) 2 



However, anticipating the result, we can replace 2/i by 
the binding energy Eb = — g— - g2 plus terms of order 
(density). We find in this wa y th at the contribution to 
ni/|A| 2 of the first term of Eq 23 after summation on k 



(m r ) 2 a = z _ x 
2tt 



(24) 



For the second term of Eq(23j due to the 5 function, the 
integration gives a factor No = miko/(2ir 2 ), the density 
of state of excess fermions at energy fi\, ti mes ( minus) 
T^(k,k-Js) at k = {k ,0}. Then we use Eq|A7] of Ap- 
pendix |A| in order to find Tg 1 up to order ko. In this way 



the contribution to 5'i 
summation on k is 



» 



|A| 2 of the second term after 



Z^NogAoll + k a AD (2tt) F aA (l,m D / mi )] (25) 



where we have introduced the dimer-particle 1 cou- 
pling strength gAD = 27T clad/itiadi where uiad 
(mimc)/(mi + ttid) is the dimer-particle 1 reduced 
mass and ciad is the dimer-fermion 1 scattering length. 
The function F a A is defined in EqA"8 of Appendix \K\ 



Concerning the third term on the r.h.s. of Eq|23| we 
use the result of Eq B2 of Appendix Ib] (k = {k, z = 
(k 2 ~k 2 )/2 mi }) 



dz 



T 3 (k, k; k) 



-Z 1 (gADm A D) 2 -rF aa 
k 



|k| m D 

fco ' TTti 



The function F a 2 is defined in Ecj B3 of Appendix [B 



After integration on k, due to the occupation number 
/(&) = ©(fco ~ |k|), we find fn <k0 $fa&n(k,k-,k). If 
we rescale k by fco we find the contribution to ni/|A| 2 
for the third term of Eq|23| 



Z- 1 g A BN koa A D rn ^ (4^ 3 ) 



mi 



X<1 



d 3 X 



7 F a 



mi 



(26) 



We just h ave to multiply the three contributions of 



Eqs |24|25 and 26 by |A| 2 in order to get the correction 



to the density of particles 1 



Sn[ a) = 



Z 



N gAD - 



IAI 



(1 + k Q a A D(Ca,l + C a ,2)) 

(27) 



The dimensionless constants C Q ,i and C a .2 are 
given by C a A = 2irF aA (l,m D /mi) and C Qj2 = 

47r 3 (m A D/TOi) J (l^i ©(! ~ x)F a , 2 (x, m D /mi). 

Going to particles of species 2, we first have to consider 
the self-energy of the particles 2 in the normal phase. 
The situation is different for £2 and Si, since, due to 
the positive chemical potential fii, the particles 2 can 
interact with the Fermi Sea of particles 1. This prob- 
lem is related to the impurity problem, which has been 
the object of various works both experimentally [13] and 
theoretically [Mj. In the normal phase, at lowest order 
in the density of particles 1, the self-energy is given by 
the so called "ladder approximation" which is shown di- 
agrammatically in Figj3j Explicitely it is 



bladder 



(k) 



Y,T2{K)G° l {K-k) 



(28) 



K 



This gives the contribution to the density of particles 2 
at lowest order 



ladder 



-E( G °( fc )) 2s 2 Qdrfer (fc) 

k 

Y,{Gl{k)) 2 T 2 (K)Gl(K-k) (29) 



k,K 



= We can make the change of variables p = K — k in 



the double summation 



-.ladder 



E k , P (G° 2 (k)) 2 T 2 (k 



p)Gl(p). In this way, we see that the integrand is ana- 
lytical for 1Ze(iu)k) < and if we perform the frequency 
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summation on iu>k, we get zero by deforming the fre- 
quency contour in the half-plane TZe(iuik) < 0. We thus 
find 



.ladder 



o 



(30) 









T3 




T3 1 



k,1 k,1 

(a) 



P.1 



k,1 



k,1 



(b) 

FIG. 2: Diagrammatic representation of s[ a \k) and s[ b \k) 



K-k,i 



k,2 



T 2 



k,2 



FIG. 3: Diagrammatic representation of Sj (ft) 



We have also the contribution of the generalized self- 
energy due to the condensate of dimers, similar to Eq|22| 
but for the particles 2 now 



Si a \k) = \A\ 2 Ti(k,k;k) 



(31) 



The contribution to the density of particles 2 is using 
Eq|I5j 



» 



J2(G° 2 (k)fsi a \k) 



(32) 



It is shown in Appendix [C] that at the order we consider 
we have (i.e. up to order (density) 5 ^ 3 ) 



» 



Z 



(33) 



Before closing this section, we notice that we could have 
written 



S 2 (k) « Z ladder (k) + S ( 2 a) (k) 
and use Eq[l2]to calculate 



(34) 



B. Gap equation 

We determine now the lowest order contribution to 
F'(k). We first have the contribution already found in 
Ref.JII] (see Figga)) 



F'( Q )(fc) 



A 



-|A|^(fc)G!^(-fc)<I>(fc,-fc;0,0) (35) 



When F'( a \k) is summed on k in the gap equation 
we recover the dimer-dimer T matrix T4 : 



20 



E 

k 



= iT 4 (0,0;0,0)|A| 2 



(36) 



However there is an other contribution to F'(k) due to 
the presence of the Fermi Sea of the extra 1 atoms, with 
one A and one hole propagator (see FigQb)) 



F'W(k) 



-G?(A)G§(-fc) £ G?(p)T 2 (p - k)Tl(k,p;p) 
v 

(37) 



Then F'^ b \k) is summed on k and the order of summa- 
tions on k and p is exchanged. The summation on k 
(including the minus sign) is precisely the integral part 
in the integral equation for T^(p,p\p) [15 . Therefore we 
have to substract the Born approximation term — Gq(— p) 



E 

k 



(%(k)G%(-k)T 2 (p-k)Tl(k,p;p) = T^p,p;p)+G° 2 (-p) 



In this way, we can transform the summation on k of 
Eq |37| in a summation on p 



E 



A 



X;gS(p)[t 3 1 (p.p;p) + g§(-p)](38) 

p 



The quantity between brackets is analytical for 
Re(ioj p ) < 0. We can therefore perform the integra- 
tion on iuj p and close the contour in the half-plane 
Re(iu! p ) < 0. We pick the contribution of the pole of 
G\(p) when |p| < fc and we get 



E 



F'( b )(fc) 



p<k, 



d 3 p 

(2t) 3 



[TZ(p,p;p) + G° 2 (-p)] 



(39) 



where p is "on the shell" p = {p, (p 2 — fe 2 )/2mi}. The 
second term in the brackets of Eq{39]is canceled by —R(0) 
in the gap equation |20| Finally, we can write a final gap 
equation in terms of the bosonic chemical potential fig 
from Eqs |20|36|39l 



1-lB 



Z\ -T 4 (0,0;0,0)|A| 2 



rrlf \ 



(40) 
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where we have omitted terms of order /j b , i.e. of the 
order of the square of density. We will see this is sufficient 
at the order we consider. 

We can express T4 in terms of the dimer-dimer 
scattering length a^D [15] . We have T 4 (0,0;0, 0) = 
Z _2 (47r/mi))(2aD£)). The dimer-dimer coupling con- 
stant is defined as usual by 9dd = (^ko-dd) /mjj. More- 
over, we can use the expansion of Appendix [A] Eq |A7| 
for T 3 1 (fc, k; k) at lowest order : T^(k, fc; k) gAo/Z. In 
this way, we can perform the integration on the wave 
vector k and we get the mean-field like contributions 
[i B = g DD (\A\ 2 /Z) + g AD (k^/6n 2 ). Using Eqsj9] and 
33 we find the bosonic chemical potential in terms of the 
densities 



Ms 



gDDn 2 + gADin x - n 2 ) 



(41) 



which is a result at linear order in density for the bosonic 
chemical potential /xg defined in Eq jl9| 



V. DENSITY EXPANSION : NEXT ORDER 

In order to get a diagram for Si at next order in fco, 
we include the possibility to have a hole line (therefore 
going backward in a graph) in addition to one A and a 
A*. Such graphs are shown in Fig{2jb). The arguments 
yielding to these graphs are given in Appendix [D] 

S[ b) (k) = \A\ 2 yTi(p,k;p)T 2 (p-k)Ti(k,p;p)Gi(p) 



(42) 



In order to calculate Snf ] = - Efc( G i( fc )) 2 - 5 i 0j ( fc )> we 
need to perform two summations on the frequencies 
and ujp. In the summation on o> p , we use the fact that 
T^(p,k;p) (— T^(k,p;p)) (this is quite easily seen by 
writing down the diagrammatic expansion where one sees 
that p always appears with a plus sign) and T 2 (p — k) 
are analytical for TZe(iuj p ) < and we perform the in- 
tegration like in Eqj23] By doing so, due to the ana- 
lytical properties of T3 1 (p, k;p) and T 2 (p — k) as a func- 
tion of the frequency iui p , we close the integration con- 
tour in lZe(iu>p) < and we get a contribution of the 
pole of G\(p) for |p| < kp and in this way we find 
that p is "on the shell" (p = {p, (p 2 - fc 2 )/2mx}). The 
summation on LUk is also easily performed if we use the 
fact that T 2 (p — k) are analytical for lZe(iuJk) > 0. 
Like in Eqj23j we get two terms due to the presence 
of a double pole in (G°(fc)) 2 , and fc is on "the shell" 
(k = {k, (k 2 — i;p)/2mi}). In the term with the deriva- 
tive of the product T^(p,k;p) 2 T 2 (k — p), we see easily 
that the derivative of T3 1 will give a finite term in the low 
density limit (as can be see easily using for instance the 
integral equation for T3 1 ). T 2 (k — p) can be written as 
the product of the bosonic propagator which is singular 
at low energy and wave vector times a regular function 
at low energy. The derivative of T 2 (k — p) at low energy 
thus gives 



\2 o( b )i 



n 2 — k 2 
2mi 



lUJ k 



(p-k) 

2m r> 



_(k2-fc2) 



D 2 -k 2 



lUJ k 



(P-k) 2 

2m r, 



(43) 



We can now perform the integrations on the wave vectors 
p and k. Like in Appendix [3] , it can be shown that the 
low wave vector domain is dominant and we get for the 
correction to the density of particles 1 



5n\ 



(*>) 



No 9 ad - 



IAI 



(k clad) (C b - C , 2 ) 



where the constant Cb is defined by 

d 3 x 0(l-|x|)(r- 



■in 



(2tt) 3 (r(l-x 2 ) + (u-x) 2 ) 



(44) 



(45) 



where u is any unitary vector. Finally, the total density 
of fermions of species 1 is given in terms of A and kg 
(= i/2toi/xi) by the expansion 



Z 

-N 



6tt 2 



IAI 



9AD- 



(1 + k a A D C(m D /m\)) 



(46) 



where the constant C{mu /mi) = C a .i +CV After calcu- 
lation of the required integrals, one finds that C is given 
in terms of the function /((m/j — mi) /(mo + mi)) given 
by Albus et al. [7]. Specifically we find it C (mo /mi) — 
(4/3)/(<$), where we have defined the dimensionless pa- 
rameter 



_ (tod - mi) 



m 2 



(m.£)+mi) 2m 1 +m 2 



(47) 



and we have 



/(^^ + i<i±fl^) M l±| )(48) 

For the particles 2, we find naturaly that there is no 
contribution of diagrams similar to the ones shown in 
Fig{2jb) (but with the roles of particles 1 and 2 ex- 
changed) because there is no hole for particles 2 at lowest 
order. The first term in the r.h.s. of Eq|46]is n 2 at the 
order we consider, due to Eq |33| Moreover we can rewrite 
Eq|46| by considering it as a Taylor expansion, after re- 
placing C by 4/(3tt)/ 



ni — n 2 



(277ii [fJ>i - 9ADn 2 (l + ^k F a A D /(#))]) 



3/2 



6tt 2 



(49) 
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FIG. 4: The two contributions to F'(fc). 



where we ignored terms of order (density) 5 I 3 . At the 
order we consider we have replaced \A\ 2 /Z by n 2 in the 
r.h.s. of the equation, and fco by Fermi wave vector kp in 
the term in factor of /. The Fermi wave vector is defined 
in terms of the density of fermions in excess in the usual 
way 



kp = (6tt (m — n 2 )) 
Eq |49| is easily inverted and we find 



1/3 



(50) 



Mi 



0A£>ra 2 I 1 H ^ — f($) (51) 



2 mi 



3tt 



and we find for the chemical potential of particules 2 



M2 



.915 15^2 



2mi 

+5a_d("-i - 2n 2 ) 

fk F a A Df(S) , 
+9 Ad I g- (3ni 



7n 2 



(54) 



where we have included terms up to order (density) 4 / 3 . 
The bosonic chemical potential ub = fJ-i + \ii + \Ej,\ is 
then from Eq[5T] and 54 



I'D 



9DD n 2 

+gAD(ni 



n 2 ) 



k F a A Df(S) 



(55) 



Looking at Eq |51| and 55 we find that we recover pre- 
cisely at this order the result of Albus et al. [7] Eqs. 
76 and 67 respectively if we consider that there are n 2 
dimers considered as pointlike bosons (mass mo) and 
Tii — n 2 fermions (mass mi). The boson-boson interac- 
tion is described by the scattering length clod and the 
fermion-bosons interaction is characterized by the scat- 
tering length a ad- 
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We can now easily calculate u 2 using Eq |51| and 
Indeed, we can follow [7] and use that Ml = 
Tp(E/V)/drii) n2 , where E is the ground state energy in 
a volume V. We integrate Eqj5l] as a function of ri\ for 
ni < Ti2, keeping n 2 fixed. We get 



E 
V 



= E B (n 2 ) 



dn' 1 ^i(n[,n 2 ) 



(52) 



where Eg(n 2 ) is a constant of integration. We see from 
this equation it is the ground state energy (per unit 
volume) of the system with an equal number of parti- 
cles 1 and 2 (" unpolarized" situation). We know from 
[IT] it is given by the Lee-Huang- Yang (LHY) [IB] ex- 
pression. Inserting Eq[5TJ into Eq 52 we get, using the 
mean field term like term of the LHY expression for 
Eb(ti 2 ) (the LHY beyond mean-field correction is of or- 
der (density) 5 ^ 2 which we neglect at the order we con- 
sider) 



E 
V 



-\E b \n 2 + -gDD(n 2 ) 
+9AD{n\ - n 2 )n 2 ( 1 



i 3 kp 
10 mi 

kpdAD 
7T 



f(S) 



(53) 



which is an expansion up to terms of order (density) 7 ^ 3 . 
From this we can take the derivative with respect to n 2 



VI. IMPURITY PROBLEM AND 
ONE-PARTICLE GAP 



For a system with N% particles 1 and N 2 particles 2, 
we denote the groundstate energy E(Nx,N 2 ). A single 
particle 2 (an "impurity") immersed in a Fermi Sea of 
particles 1 has an excitation energy 

E(Ni,l)-E(Ni,0) = M2(ni,n 2 = 0) (56) 

In the BEC limit, this problem has been the subject of 
many studies recently both experimentally |13) and the- 
oretically [13] • From Eq 54 with n 2 = 0, we find 



vmp 

M 2 



kp 
1tri\ 



\E h 



+9 ad ni 1 



kp CLADf(d) 
7T 



(57) 



where kp = (67r 2 ni) 1 / 3 . This result, and also a next- 
order correction, was found by Mora and Chevy in 
Ref.pU. 

The one-particle gap A ffop is related to the ground- 
state energy via [10], [17] 

A gap = ~ (2 E(N + 1,N) — E(N + 1, iV + 1) — E(N, N)) 



(Mi - M2) 



(58) 



From Eq 51|54[ we find 

Agap = ^(\E b \ + (2g AD -g DD )n 2 ) (59) 

a result given in Eq.(40) (for particles of equal masses) 
ofRef.QZJ]. 



VII. COMMENTS AND CONCLUSION 



from T3 1 by setting \i\ to and ^2 to — \E\\. One can 
define in the same way T 2v , G\ and G\ . We write 
Ti(k, k; k) - Tl v (0, 0; 0) = (T^k, k; k) - k; k j) + 

{Tl v (k, k; k) - 7^(0, 0; 0)), and we calculate the two dif- 
ferences between parentheses separately. 



Ti(k,k;k)-Tl v (k,k;k) 



We have presented in this work a low density expan- 
sion of the equation of state of a condensate of compos- 
ite bosons made of two fermions of two different species, 
in the case where one of the two species is in excess. 
We have shown results of the expansion for the chemi- 
cal potentials up to order (density) 4 ^ 3 . However there 
must be terms of the type of the LHY correction of order 
(density) 3 ^ 2 (more precisely of the type v^/ 2 and most 

1 /2 

probably of type nin 2 ' ). In order to see the quantita- 
tive importance of the LHY correction and the correction 
to [1b given in Eq[55j we calculate the ratio R between 
these two terms 



R = 



S2 / ^ \ 1/2 

9AD(ni - n 2 )k F a AD typ- 



0.258 



"2 

np 



4/3 



F,b a 



(60) 



with uf = n\ — ri2 and kp.b = (6n 2 n2) 1 ^ 3 . We have 
taken m\ = m 2 , ap>p)/a = 0.6 [THUS], o-ad/ - = 1-18 
|19j in order to get the numerical prefactor. In principle, 
for fermions densities of the same order of magnitude, 
the ratio 712/ 'uf is of order 1, and the ratio R tends to 
in the dilute limite ^kpi, a 0. However, in the 
experiment of Ref.jl, looking at Fig. 4a, we estimate at 
the center of the trap TI2I tif ~ 8 and kp^a « 1. This 
gives R ~ 4.1. Close to the edge of the bosonic cloud, if 
we estimate that the dimer density is divided by 2, this 
gives ^Jkp^a « 0.9 and for a ratio n 2 /riF ~ 4, we find 
R w 1.46. From this order of magnitude estimate, we see 
that including the corrections of Eqs {5~l|54| and[55"|for the 
equations of state (in addition to the LHY correction as 
it was done in Ref.pQ) should improve quantitatively the 
analysis of the experiment of Ref. pQ . 
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Appendix A: Corrections to T3 1 

We want first to determine at lowest order in the 
densities the difference T^(k,k;k) — T^ v (0, 0;0) where 
k = {k,(k 2 - fc2)/(2mi)} and |k| < a~\ T^ v is 
the dimer-atom 1 T-matrix in vacuum. It is obtained 




-qi,2 k-qi p k-q p ^ k-q N -q N ,2 




qi,1 ™ q P ,1 ™ qN,1 



FIG. 5: The diagrammatic expansion of T 3 (k, k; k). 



In order to calculate the first difference, we note that it 
comes from the fact that the fermionic propagators and 
the dimer propagator are different in vacuum and with 
Hi 7^ and 1x2 7^ — | | . The principle of the method is 
therefore the following : we write AGi^fc) = G\ 2 (k) — 
G v ia {k), AT 2 {K) = T 2 (K) - T 2 , V {K) and we expand the 
diagrammatic series for T3 1 in powers of AGi^ and AT 2 - 
If we look at the diagrams for T 3 1 (fc, k; k) with N dimer 
propagators T 2 , we see that G° appear N times in G\(q p ) 
forp = 1,- • -N (see Figj|). If we write G%(q p ) = G v 1 (q p ) + 
AGi(q p ), the term of order AG\(q p ) will be multiplied 
by T 2tV (k — q p ), by all the diagrams of v (k,k;k) on 
his left and by all the diagrams of T3 1 v (k,k;k) on his 
right. However those diagrams are Feynman diagrams 
for T^ v (k, q; k) and T^ v (q,k;k) respectively. Therefore 
when we resum all the diagrams with any number of T 2 
we have 

Y, T L(k, q; k)T 2 , v (k - q)AG 1 (q)Tl v (q, k; fc)(Al) 



where we have taken the atom-atom and atom-dimer T 
matrices in vacuum at lowest order. By looking at the 
diagrammatic expansion of T3 v (k, q; k) (= T$ v (q, k;k)), 
we find that it is analytic for TZe(iu> q ) > (q = {q, iui q } 
appears always with a minus sign in a particle 2 prop- 
agator G2). For the same reason T2(k — q) is analytic 
for 7te(iu)q) > 0. In the zero temperature limit, we can 
therefore transform the frequency summation in Eq Al 
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in an integral along the imaginary axis. Then we deform 
this contour in the half-plane 1Ze(iu) q ) > 0. We get the 
two contributions of the poles of AG\ (q) at the frequen- 
cies (q 2 -£; 2 )/2mi and q 2 /2mi ("on the shell") for G\(q) 
and Gl(q) respectively 



^[(Tl v (k,q-,k)) 2 T 2iV (k-q)e(\q\-k ) 
-(Tl v (k,q;k)) 2 T 2 , v (k-q))} 



(A2) 



where q = {q, (q 2 — fc 2 )/2mi} and q = {q, q 2 /2mi} for 
the first and second terms respectively. The main con- 
tribution to the integral comes from small wavevectors 
|q| <C a -1 as is explained below. In this case, we can re- 
place Ts <v by its low energy value Ts <v {0, 0; 0) = Z~ 1 gAD 
(the Z factor comes from the two external dimer lines, 
each one giving a \/~Z factor [T51 120) ). For such small wave 
vectors, we can also take T 2<V ({IC, iQkY) ~ Z/(iflx — 
K 2 /2m£>) (a bosonic propagator multiplied by two V~Z 
factors for the incoming and outgoing dimer end). We 
therefore get at lowest order 

Ti(k, k; k) - Th, (k, k; k) = Z-\g AD fy. 



1 



d 3 q 

(2tt) 3 



e(|q| - fc ) 
(k 2 -q 2 ) _ (k-q) 



2m! 



2m r 



(k^-fcg-q^) _ (k-q)2 



2m i 



2m i 



(A3) 



We note that this correction is of order ko (since |k| is 
of order ko), i.e. of order (density) 1 / 3 , due to the "small 
denominators" of order k% . 

We show now that the contribution to integral A2 
for wavevectors |q| > A, where A is a cut-off, with 
to « A « a -1 is negligible. Indeed, in such a range of 
wavevectors, we can perform a Taylor expansion of the 



difference in Eq A2 with respect with the energy variable 
of q. The energy difference is of order fc 2 , and by close 
inspection of the high wavectors behavior of and T 2 , 
we find a converging integral at infinity. The low wave 
vector behavior (~ A) is dominated by the small denom- 
inators of the T 2 's and this gives a contribution of order 



&q A. Compare to the term of Eq A3 this is negligible if 



koAa 2 <C 1. This is possible if A = (fco/a) 1 / 2 for instance. 
As a conclusion, the high wavevector domain gives a neg- 
ligible contribution. 

We can make the same sort of analysis with AT 2 . In 
this case, we would find that the characteristic low wave 
vector is of order the square root of the density (this 
is obtained by comparing \ig and R to K 2 in the de- 
nominator of Eql6). This gives a contribution of the 
order of (density) 1 ' rz instead of ko which is of the order 
of (density) 1 ^ 3 . Therefore it is negligible. 

We can also consider the effect of AG 2 . There is no 
"small denominator" in this case, since one has to com- 
pare the kinetic energy q 2 /2m 2 to the (large) chemical 
potential [i 2 ~ — \Eb\. As a result, we find the correction 
is of order the chemical potential correction n 2 + \Ef,\, 
which is of order fc 2 . Therefore it is negligible here. 



2. T 3 1 ,„(fc,fc;fc)-r 3 1 1 ,(0,0;0) 

In order to calculate this difference, wc note that in the 
diagrammatic expansion of T 3 1 , u (fc, k; k), the 4- wavevector 
k appears only in G 2 (k — q p — q p +i) and in T 2 (k — k p ) 
(not in the Gfs). We define now AT 2 (k - q p ) = T 2 (k - 
q p ) — T 2 (—q p ) and similarly for AG 2 . 

The corrections of order AG 2 gives a negligible term 
again, for the same reason as given in subsection | A I| 

The correction of order AT 2 can be analyzed like the 
correction of order AG\ above. Indeed, we can permute 
the roles of T 2 (k — q p ) and G\(q p ) in the argument. In 
this way, we find 



r 3 y0, q; Q)AT 2tV (k - q)G v 1 (q)Tl v (q, 0; 0)(A4) 



The frequency summation is done as above and we get 
the contribution "on the shell" 



d3f l ,ml r_ ^ 1 . n\\2 It\ /n . _ k 2 - fc^ - q 2 



(27T) 



(T 3 %(0,{q,^};0)) 2 (T 2 ^({k-q,. 



2toi 



2m! 



}) 



-T 2)t ,({-q,- 



2m 



■}) 



(A5) 



We can show again that the low wave vector domain gives 
the main contribution, and replace the atom-dimer T ma- 
trix by Z~ x gAD and the T 2 's by their bosonic expressions. 
In this way we get at lowest order 



TUk,k;k) - Z-^ad = Z- l (g AD f 



d 3 q 



(2^) E 



(k 2_ fc 2_ q2) 

2m i 



(k-q) 2 
2m n 



2m A L 



(A6) 



Collecting Eq\M\ and \M\ we get 



T^(k,k;k) - Z 1 g AD = 
Z- l (g AD ) 2 (k )m AD F aA (g, aa) (A7) 

where the dimensionless function F aj i is defined by 



F a A^r) = 2/1*3 



e(|y|-l)(r + l) 
r ( x 2 _ y 2) _ ( x - y) 2 



1 

V 2 

(A8) 



Appendix B: Evaluation of ~Tl(k,k;k) 

In order to evaluate this derivative, we note again that 
k appears in the diagrammatic expansion of T 3 1 (fc,fc;/c) 



10 



in the propagators of the particles 2 G 2 (k — q p — q p +i) 
and in the dimer propagators T 2 (k — q p ). If we take the 
derivative of T 2 (k — q) with respect to the frequency of 
k, which we denote z here, we get similarly to |A 2| 



Ti(k,q;k)—T 2 (k 



q)G° 1 (q)Ti{q,k;k) (Bl) 



We can perform the frequency summation as before, 
which has the effect to take the frequency "on the shell" 
(iuj q = (q 2 — fc 2 )/2mi) and |q| > k . Again, the 
low wave-vector domain gives the main contribution. 
We can safely replaces the Tg's by Z~ x qad and T 2 by 
the bosonic expression. We get at lowest order (using 
k = {k, (k 2 - fc 2 )/2m!}) 



d 

dz 



Tl{k,k;k) = -Z 1 (gAD'm AD ) 2 



-i^f — ,— ' 
ko a ' 2 V &o ' mi 

(B2) 



where 



F a , 2 {x,r) = 4 



ri 3 y 8(|y|-l)(r + l) 2 
(2tt)3 ( r (x 2 - y 2 ) - (x - y) 2 f 



(B3) 



Appendix C: Calculation of n 



(a) 



We proceed like for the calculation of correction to T3 : 
for each propagator (G?, G 2 or T 2 ) we write it as the 
sum of its value in vacuum plus a difference. In this 
way, we are left with most of the propagators in vacuum. 
Therefore we start with all diagrams in vacuum for ■ 
The calculation of the contribution to n 2 is then similar 
to Ref. [TT] , due to the analytical properties of T 2 V (k, k; k) 

I Al 2 

as a function of the frequency. We find n 2 = ■ In the 
course of the calculation, we have used 2fi s» —\E^\ and 
we have neglected terms of the order of the density. 
We turn now to corrections to propagators. 



Correction of order AGi 



We first consider corrections of order AGi . Looking 



at the diagrams for they look like Fig|5| but with 

the roles of particles 1 and 2 exchanged. We see that 
Gi appears in four places : i) G\(—q\), ii) G\{— q^), in) 
G\(k — q p — q p +i) and iv) G\(—k). 

In case i), we find the contribution to n 2 

|A| 2 ^(G l 2 (fc)) 2 AG 1 (- <Zl )T 2 (A : - qi )GUqi)Tl v (q 1} k; k) 

k,qi 

But in the diagrammatic expansion of T3 v (qi, k; k), the 
variable k only appears with a plus sign. This means 
T^ v (q 1 ,k;k) is analytical for lZe(ioJk) < 0. This is also 



true for T 2 (k — qi), therefore, when we sum on ujk, we 
find zero. The case ii) is similar if we replace q± by qjy, 
and we also find a result equal to zero. 

The case Hi) is different. We find, similarly to what was 
found in |A 1[ the following contribution to the density of 
particles 2 

E (G 2 (k)) 2 Tl v (k, q; k)T 2 (k q)Gl{q) 

k,q,q' 

AGt(k - q - q')Tl v (q', k; k)T 2 (k - q')Gl{q') (CI) 

All the integrand except AGi(fc — q — q') is analytical 
for 1Ze{iuJk) < 0. As usual by deforming the contour of 
integration in the half plane 1Ze{iuJk) < 0, we get a pole 
of G\{k — q — q') for iwk = iu q + iuj q > + ((k — q — q') 2 — 
fc 2 )/(2mi) if |k - q - q'| < k . The difference ((k - q - 
q') 2 — fcg)/ (2toi) is small compared to the binding energy 
\Eb |, which is the energy scale of T^ v . Therefore we can 
replace at lowest order iuik by iuj q + iuj q > in Eq CI In the 



same manner, since |k — q— q' 



< fc is much smaller than 
a " a typical wave vector scale for T^ v , we can replace 
k by q + q' in Eq CI This means that we can replace 
everywhere, at lowest order, the 4-momentum k by q + q' 
in Eq |Cl| The integrand is then independant of k and 
we get a phase space factor /cq/(67t 2 ). The remaining 
summations on q and q' give the factor 

E( G 2 (9 + q')) 2 TL(q + q', q; q + q')T 2 (q')G v 2 ( q ) 
9,3' 

TIM, q + q';q + q')T 2 {q)G v 2 ( q ') (C2) 

We can perform the summation on the frequencies as 
before (only the "Born" contributions to Tf v contribute) 
and we get the two contributions "on the shell" iuj q = 
— q 2 /(2mi) and iu! q > = — q' 2 /(2mi) 

GS(q', -~) I g v M + q', -^^) ) (C3) 



2mi 



2toi 



This integral is convergent close to q = and q' = 0. It is 
also easily shown that it converges when the wave vectors 
tend to infinity. The integral [C3| is therefore a constant, 
independant of k . As a conclusion, this case Hi) gives 
a contribution to n 2 of order |A| 2 (fco) 3 = (density) 2 and 
we can safely neglect it. 

Finally we consider the case iv). This is given by 



|A| 2 E(^(fc)) 2 AG 1 (-fc) 



(C4) 



The frequency integration is easily done and we find 



e(fc-fcp) 



(2ir) 3 



(C5) 
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Performing the integral on k, we find the contribution to 
n 2 

m\A\ 2 Tr(m r a) 3 (C6) 
2. Correction of order AG2 



Conclusion 



As a conclusi on of th is Appendix, we find that the two 
corrections Eqs C6|C9 of order |A| 2 £;q exactly compen- 
sate, and we find n 2 — \A\ 2 /Z and there is no correction 
up to order (density) 5 ^ 3 included. 



Next, we study the contributions of order AG2. In a 
diagram for (G^fc)) 2 ^ (fc), a particle 2 propagator can 
be i) in the two external propagators G°(fc) or «,) in the 
dimer-(fermion 2) T-matrix T 2 (k,k;k). In case ij the 
correction to n 2 of order AG2 is given by 

- 2|A| 2 J2 G v 2 (k)AG 2 (k)Tl v (k, fc; fc) (C7) 
k 



We can then write T^ v as the sum of 



G\(-k), which 
has a pole for — k 2 /(2mi) < 0, plus a function which is 
analytical for lZe(iuik) < [IT]. The frequency summa- 
tion is then done as usual by deforming the contour of 
inegration in the half plane lZe(iujk) < 0. We get in this 
way, keeping terms of order (density) 2 ^ 3 |A| 2 



\E h 



d 3 k 



k 2 +q- 
2m r 



(C8) 



But /Z2 + l-Ebl = /is — Mi = ^ Mi + O(density). And at 
lowest order we can replace — /i2 by . This gives the 
contribution 

-Aii|A| 2 7r(m r a) 3 (C9) 

The case n) is similar to what we discussed in [XT] ex- 
cept that the roles of particles 1 and 2 are exchanged. 
Therefore we find for the contribution to n 2 



(Gm) 2 T 2 (k, q; k)T 2 (k - q)AG 2 (q)T 2 Jq, k; k) 



k , q 



(CIO) 



But all the functions are analytical for 1Ze(iujk) < 0. 
Therefore, we can perform the summation on u)^ by 
transforming it into an integral, and deform the contour 
in the half-plane TZe(icJk) < 0. In this way, we find that 
the sum of Eq |C10| is zero. 



3. Correction of order AT2 

We finish this appendix by adressing the contributions 
to n 2 of order AT2. They are similar to the ones of 
order AG2, since for every dimer propagator T 2 there is 
a fermion 2 propagator G2. We find 



(G v 2 (k)rTl v (k, q; k)AT 2 (k - q)G v 2 (q)Tl v (q, k; k) 



k . q 



(Cll) 

Similarly to Eq jC10[ we find that this contribution is 



Appendix D: Feynman diagrams of S{ 




FIG. 6: (a) A general diagram for Si(fc) with one hole line, 
one A and one A*, (b) Diagrammatic expansion of the new 
vertex. Only the first diagram contributes at lowest order (see 
text). 



We explain here why the Feynman diagrams for the 
generalized self-energy S%(k) are given by Figgb) at the 
order we consider. We look for graphs with two anoma- 
lous self-energies (one A, one A*) and one hole line. In 
general this can be drawn as shown in Fig[6](a). The 
white box is a new vertex with one incoming dimer (made 
of a fermion of species 1 and a fermion of species 2) , two 
incoming fermions of species 1 and one outgoing dimer 
and two outgoing fermions of species 1. This sort of ver- 
tex is very similar to the vertex $ introduced in [15] 
with one incoming dimer and two incoming fermions, 
which do not interact first (this condition is automati- 
cally fulfilled here since the two incoming fermions are 
of the same species). In this case, one of the incoming 
fermions 1 first interacts with the fermion 2 of the incom- 
ing dimer. All the subsequent processes are resummed in 
the dimer-fermion T-matrix T 3 X and we have one dimer 
and two fermions 1 at the output of T 3 X . Then the other 
fermion 1 interacts with the dimer, giving an other T3 1 , 



(see Fig 
in Fig|6 



6[b)) etc ... Consider the first term with two T 3 s 



b). In order to get the generalized self-energy, 
we have to close a fermion 1 loop. It is then rather easy to 
see that in order to avoid double counting, we must con- 
sider the diagram of Fig 2jb). As is shown in Eq 44 this 
gives a contribution to the density of fermions 1 of order 
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|A| 2 fc§. This order of magnitude can be estimated easily 
: the condition |k| < k gives a factor fcg, the hole line in 
Si gives an other factor k$. But there is a large term of 
order 1/fcg coming from the derivative of T^p — k). This 
is so because both p and k are at low energy and therefore 
there is a small energy denominator of order feg/(2mi). If 



we consider now the other terms with three or more Tg 1 , 
we will get the same phase space factor (ko) 6 , but as it is 
easily seen by writing the diagrams, the dimer propaga- 
tors will never be at low energy, and we will not get the 
large prefactor l/(fcg) 4 . This concludes the justification 
of retaining diagrams of Fig[2jb). 
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